Glycogen phase separation drives macromolecular rearrangement and asymmetric division in E. coli.

Bacteria often experience nutrient limitation in nature and the laboratory. While exponential and stationary growth phases are well characterized in the model bacterium Escherichia coli , little is known about what transpires inside individual cells during the transition between these two phases. Through quantitative cell imaging, we found that the position of nucleoids and cell division sites becomes increasingly asymmetric during transition phase. These asymmetries were coupled with spatial reorganization of proteins, ribosomes, and RNAs to nucleoid-centric localizations. Results from live-cell imaging experiments, complemented with genetic and 13 C whole-cell nuclear magnetic resonance spectroscopy studies, show that preferential accumulation of the storage polymer glycogen at the old cell pole leads to the observed rearrangements and asymmetric divisions. In vitro experiments suggest that these phenotypes are likely due to the propensity of glycogen to phase separate in crowded environments, as glycogen condensates exclude fluorescent proteins under physiological crowding conditions. Glycogen-associated differences in cell sizes between strains and future daughter cells suggest that glycogen phase separation allows cells to store large glucose reserves without counting them as cytoplasmic space.

Epidemiology and genetic diversity of invasive Neisseria meningitidis strains circulating in Portugal from 2003 to 2020.

Invasive meningococcal disease (IMD) continues to be a public health problem due to its epidemic potential, affecting mostly children. We aimed to present a detailed description of the epidemiology of IMD in Portugal, including insights into the genetic diversity of Neisseria meningitidis strains. Epidemiological analysis included data from the Portuguese National Reference Laboratory of Neisseria meningitidis during 2003 to 2020. Since 2012, N. meningitidis isolates have also been assessed for their susceptibility to antibiotics and were characterized by whole genome sequencing. During 2003-2020, 1392 confirmed cases of IMD were analyzed. A decrease in the annual incidence rate was observed, ranging from 1.99 (2003) to 0.39 (2020), with an average case fatality rate of 7.1%. Serogroup B was the most frequent (69.7%), followed by serogroups C (9.7%), Y (5.7%), and W (2.6%). Genomic characterization of 329 isolates identified 20 clonal complexes (cc), with the most prevalent belonging to serogroup B cc41/44 (26.3%) and cc213 (16.3%). Isolates belonging to cc11 were predominantly from serogroups W (77.3%) and C (76.5%), whereas cc23 was dominant from serogroup Y (65.7%). Over the past 4 years (2017-2020), we observed an increasing trend of cases assigned to cc213, cc32, and cc11. Regarding antimicrobial susceptibility, all isolates were susceptible to ceftriaxone and 61.8% were penicillin-nonsusceptible, whereas 1.4% and 1.0% were resistant to ciprofloxacin and rifampicin. This is the first detailed study on the epidemiology and genomics of invasive N. meningitidis infections in Portugal, providing relevant data to public health policy makers for a more effective control of this disease.

Factors Associated with the Extent of Clinical Attachment Loss in Periodontitis: A Multicenter Cross-Sectional Study.

Periodontitis has significant public health implications, affecting individuals' overall health, well-being, and quality of life. This study aimed to assess the risk factors associated with the extent of clinical attachment loss (CAL) in a population diagnosed with periodontitis. Six hundred and sixty-seven patients with different degrees of CAL (mild, n = 223; moderate, n = 256; and advanced, n = 188) were enrolled. Socio-demographics, lifestyle, microbiological profiles, specific immune response, obesity, and single-nucleotide polymorphism of the IL1 gene were determined. Unconditional logistic regression models were conducted to determine the factors associated with the extent of CAL. Aging, smoking, microbial factors, plaque index, and IgG2 antibodies against Aggregatibacter actinomycetemcomitans were associated with advanced CAL. IgG2 antibodies against A. actinomycetemcomitans (OR 1.50; CI 95% 1.23-1.81), plaque accumulation (OR 2.69; CI 95% 2.20-3.29), Porphyromonas gingivalis (OR 1.93; CI 95% 1.35-2.76), Tanerella forsythia (OR 1.88; CI 95%1.30-2.70), and current smoking (OR 1.94; CI 95% 1.31-2.87) were associated with advanced CAL. Gene IL polymorphisms, obesity, and stress were not associated with the extent of CAL. Aging, plaque accumulation, smoking, and having antibodies against A. actinomycetemcomitans were the most critical factors associated with advanced CAL. In contrast, obesity, stress, and gene polymorphisms were not associated with the extent of CAL.

Genome-Scale Characterization of Mycobacterium abscessus Complex Isolates from Portugal.

The Mycobacterium abscessus complex (MABC) is an emerging, difficult to treat, multidrug-resistant nontuberculous mycobacteria responsible for a wide spectrum of infections and associated with an increasing number of cases worldwide. Dominant circulating clones (DCCs) of MABC have been genetically identified as groups of strains associated with higher prevalence, higher levels of antimicrobial resistance, and worse clinical outcomes. To date, little is known about the genomic characteristics of MABC species circulating in Portugal. Here, we examined the genetic diversity and antimicrobial resistance profiles of 30 MABC strains isolated between 2014 and 2022 in Portugal. The genetic diversity of circulating MABC strains was assessed through a gene-by-gene approach (wgMLST), allowing their subspecies differentiation and the classification of isolates into DCCs. Antimicrobial resistance profiles were defined using phenotypic, molecular, and genomic approaches. The majority of isolates were resistant to at least two antimicrobials, although a poor correlation between phenotype and genotype data was observed. Portuguese genomes were highly diverse, and data suggest the existence of MABC lineages with potential international circulation or cross-border transmission. This study highlights the genetic diversity and antimicrobial resistance profile of circulating MABC isolates in Portugal while representing the first step towards the implementation of a genomic-based surveillance system for MABC at the Portuguese NIH.

Viral genetic clustering and transmission dynamics of the 2022 mpox outbreak in Portugal.

Pathogen genome sequencing during epidemics enhances our ability to identify and understand suspected clusters and investigate their relationships. Here, we combine genomic and epidemiological data of the 2022 mpox outbreak to better understand early viral spread, diversification and transmission dynamics. By sequencing 52% of the confirmed cases in Portugal, we identified the mpox virus sublineages with the highest impact on case numbers and fitted them into a global context, finding evidence that several international sublineages probably emerged or spread early in Portugal. We estimated a 62% infection reporting rate and that 1.3% of the population of men who have sex with men in Portugal were infected. We infer the critical role played by sexual networks and superspreader gatherings, such as sauna attendance, in the dissemination of mpox virus. Overall, our findings highlight genomic epidemiology as a tool for the real-time monitoring and control of mpox epidemics, and can guide future vaccine policy in a highly susceptible population.

Candida auris in Intensive Care Setting: The First Case Reported in Portugal.

Candida auris is an opportunistic human pathogen that has rapidly spread to multiple countries and continents and has been associated with a high number of nosocomial outbreaks. Herein, we report the first case of C. auris in Portugal, which was associated with a patient transferred from Angola to an ICU in Portugal for liver transplantation after a SARS-CoV-2 infection. C. auris was isolated during the course of bronchoalveolar lavage, and it was subjected to antifungal susceptibility testing and whole-genome sequence analysis. This isolate presents low susceptibility to azoles and belongs to the genetic clade III with a phylogenetic placement close to African isolates. Although clade III has already been reported in Europe, taking into account the patient's clinical history, we cannot discard the possibility that the patient's colonization/infection occurred in Angola, prior to admission in the Portuguese hospital. Considering that C. auris is a fungal pathogen referenced by WHO as a critical priority, this case reinforces the need for continuous surveillance in a hospital setting.

Development of an amplicon-based sequencing approach in response to the global emergence of mpox.

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Comparative complete scheme and booster effectiveness of COVID-19 vaccines in preventing SARS-CoV-2 infections with SARS-CoV-2 Omicron (BA.1) and Delta (B.1.617.2) variants: A case-case study based on electronic health records.

Background: Information on vaccine effectiveness in a context of novel variants of concern (VOC) emergence is of key importance to inform public health policies. This study aimed to estimate a measure of comparative vaccine effectiveness between Omicron (BA.1) and Delta (B.1.617.2 and sub-lineages) VOC according to vaccination exposure (primary or booster). Methods: We developed a case-case study using data on RT-PCR SARS-CoV-2-positive cases notified in Portugal during Weeks 49-51, 2021. To obtain measure of comparative vaccine effectiveness, we compared the odds of vaccination in Omicron cases versus Delta using logistic regression adjusted for age group, sex, region, week of diagnosis, and laboratory of origin. Results: Higher odds of vaccination were observed in cases infected by Omicron VOC compared with Delta VOC cases for both complete primary vaccination (odds ratio [OR] = 2.1; 95% confidence interval [CI]: 1.8 to 2.4) and booster dose (OR = 5.2; 95% CI: 3.1 to 8.8), equivalent to reduction of vaccine effectiveness from 44.7% and 92.8%, observed against infection with Delta, to -6.0% (95% CI: 29.2% to 12.7%) and 62.7% (95% CI: 35.7% to 77.9%), observed against infection with Omicron, for complete primary vaccination and booster dose, respectively. Conclusion: Consistent reduction in vaccine-induced protection against infection with Omicron was observed. Complete primary vaccination may not be protective against SARS-CoV-2 infection in regions where Omicron variant is dominant.

Comparative Effectiveness of COVID-19 Vaccines in Preventing Infections and Disease Progression from SARS-CoV-2 Omicron BA.5 and BA.2, Portugal.

We estimated comparative primary and booster vaccine effectiveness (VE) of SARS-CoV-2 Omicron BA.5 and BA.2 lineages against infection and disease progression. During April-June 2022, we implemented a case-case and cohort study and classified lineages using whole-genome sequencing or spike gene target failure. For the case-case study, we estimated the adjusted odds ratios (aORs) of vaccination using a logistic regression. For the cohort study, we estimated VE against disease progression using a penalized logistic regression. We observed no reduced VE for primary (aOR 1.07 [95% CI 0.93-1.23]) or booster (aOR 0.96 [95% CI 0.84-1.09]) vaccination against BA.5 infection. Among BA.5 case-patients, booster VE against progression to hospitalization was lower than that among BA.2 case-patients (VE 77% [95% CI 49%-90%] vs. VE 93% [95% CI 86%-97%]). Although booster vaccination is less effective against BA.5 than against BA.2, it offers substantial protection against progression from BA.5 infection to severe disease.

Molecular Capture of Mycobacterium tuberculosis Genomes Directly from Clinical Samples: A Potential Backup Approach for Epidemiological and Drug Susceptibility Inferences.

The application of whole genome sequencing of Mycobacterium tuberculosis directly on clinical samples has been investigated as a means to avoid the time-consuming need for culture isolation that can lead to a potential prolonged suboptimal antibiotic treatment. We aimed to provide a proof-of-concept regarding the application of the molecular capture of M. tuberculosis genomes directly from positive sputum samples as an approach for epidemiological and drug susceptibility predictions. Smear-positive sputum samples (n = 100) were subjected to the SureSelectXT HS Target Enrichment protocol (Agilent Technologies, Santa Clara, CA, USA) and whole-genome sequencing analysis. A higher number of reads on target were obtained for higher smear grades samples (i.e., 3+ followed by 2+). Moreover, 37 out of 100 samples showed ≥90% of the reference genome covered with at least 10-fold depth of coverage (27, 9, and 1 samples were 3+, 2+, and 1+, respectively). Regarding drug-resistance/susceptibility prediction, for 42 samples, ≥90% of the >9000 hits that are surveyed by TB-profiler were detected. Our results demonstrated that M. tuberculosis genome capture and sequencing directly from clinical samples constitute a potential valid backup approach for phylogenetic inferences and resistance prediction, essentially in settings when culture is not routinely performed or for samples that fail to grow.

Development of an amplicon-based sequencing approach in response to the global emergence of human monkeypox virus.

The 2022 multi-country monkeypox (mpox) outbreak concurrent with the ongoing COVID-19 pandemic has further highlighted the need for genomic surveillance and rapid pathogen whole genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for SARS-CoV-2. Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical samples that tested presumptive positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR cycle threshold below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples.

Introduction: Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited. Methods: To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity. Results: We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries. Discussion: Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.

Global gene expression analysis of Streptococcus agalactiae at exponential growth phase.

A comparative overview of the global gene expression levels of S. agalactiae reference strain NEM316 at the exponential growth phase was done through RNA-sequencing. The expression levels of 47 genes potentially linked to virulence evidenced that: i) the major nuclease, GBS_RS03720/gbs0661, presented higher mean expression values than the remainder of DNase genes; ii) the genetic pilus island PI-2a genes presented higher mean expression values than PI-1 coding genes; and, iii) three virulence-associated genes ranked among the top-100 most expressed genes (GBS_RS07760, GBS_RS09445 and GBS_RS03485). Among this top-100, genes encoding proteins involved in "Translation, ribosomal structure and biogenesis" represented 46%. Curiously, genes with no assigned function were grouped in the category of highly expressed genes. As very little is known about the molecular mechanisms behind the release of DNases, preliminary assays were developed to understand whether direct DNA exposure would affect gene expression at the exponential growth phase. No differentially expressed genes were detected, indicating that follow-up studies are needed to disclose the complex molecular pathways (and stimuli) triggering the release of DNases. In general, our insights on the global expression levels of NEM316 at exponential growth phase with and without DNA exposure should open novel research lines to decipher S. agalactiae puzzling adaptation and virulence mechanisms, such as DNase production.

Bioinformatics study of expression from genomes of epidemiologically related MRSA CC398 isolates from human and wild animal samples.

One of the most important livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) genetic lineages is the clonal complex (CC) 398, which can cause typical S. aureus-associated infections in people. In this work, whole-genome sequencing, RNA-sequencing, and gel-based comparative proteomics were applied to study the genetic characteristics of three MRSA CC398 isolates recovered from humans (strains C5621 and C9017), and from an animal (strain OR418). Of the three strains, C9017 presented the broadest resistance genotype, including resistance to fluroquinolone, clindamycin, tiamulin, macrolide and aminoglycoside antimicrobial classes. The scn, sak, and chp genes of the immune evasion cluster system were solely detected in OR418. Pangenome analysis showed a total of 288 strain-specific genes, most of which are hypothetical or phage-related proteins. OR418 had the most pronounced genetic differences. RNAIII (δ-hemolysin) gene was clearly the most expressed gene in OR418 and C5621, but it was not detected in C9017. Significant differences in the proteome profiles were found between strains. For example, the immunoglobulin-binding protein Sbi was more abundant in OR418. Considering that Sbi is a multifunctional immune evasion factor in S. aureus, the results point to OR418 strain having high zoonotic potential. Overall, multiomics biomarker signatures can assume an important role to advance precision medicine in the years to come. SIGNIFICANCE: MRSA is one of the most representative drug-resistant pathogens and its dissemination is increasing due to MRSA capability of establishing new reservoirs. LA-MRSA is considered an emerging problem worldwide and CC398 is one of the most important genetic lineages. In this study, three MRSA CC398 isolates recovered from humans and from a wild animal were analyzed through whole-genome sequencing, RNA-sequencing, and gel-based comparative proteomics in order to gather systems-wide omics data and better understand the genetic characteristics of this lineage to identify distinctive markers and genomic features of relevance to public health.

Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus.

The largest monkeypox virus (MPXV) outbreak described so far in non-endemic countries was identified in May 2022 (refs. 1-6). In this study, shotgun metagenomics allowed the rapid reconstruction and phylogenomic characterization of the first MPXV outbreak genome sequences, showing that this MPXV belongs to clade 3 and that the outbreak most likely has a single origin. Although 2022 MPXV (lineage B.1) clustered with 2018-2019 cases linked to an endemic country, it segregates in a divergent phylogenetic branch, likely reflecting continuous accelerated evolution. An in-depth mutational analysis suggests the action of host APOBEC3 in viral evolution as well as signs of potential MPXV human adaptation in ongoing microevolution. Our findings also indicate that genome sequencing may provide resolution to track the spread and transmission of this presumably slow-evolving double-stranded DNA virus.

SARS-CoV-2 introductions and early dynamics of the epidemic in Portugal.

Background: Genomic surveillance of SARS-CoV-2 in Portugal was rapidly implemented by the National Institute of Health in the early stages of the COVID-19 epidemic, in collaboration with more than 50 laboratories distributed nationwide. Methods: By applying recent phylodynamic models that allow integration of individual-based travel history, we reconstructed and characterized the spatio-temporal dynamics of SARS-CoV-2 introductions and early dissemination in Portugal. Results: We detected at least 277 independent SARS-CoV-2 introductions, mostly from European countries (namely the United Kingdom, Spain, France, Italy, and Switzerland), which were consistent with the countries with the highest connectivity with Portugal. Although most introductions were estimated to have occurred during early March 2020, it is likely that SARS-CoV-2 was silently circulating in Portugal throughout February, before the first cases were confirmed. Conclusions: Here we conclude that the earlier implementation of measures could have minimized the number of introductions and subsequent virus expansion in Portugal. This study lays the foundation for genomic epidemiology of SARS-CoV-2 in Portugal, and highlights the need for systematic and geographically-representative genomic surveillance.

Unraveling the hurdles of a large COVID-19 epidemiological investigation by viral genomics.

COVID-19 local outbreak response relies on subjective information to reconstruct transmission chains. We assessed the concordance between epidemiologically linked cases and viral genetic profiles, in the Baixo Vouga Region (Portugal), from March to June 2020. A total of 1925 COVID-19 cases were identified, with 1143 being assigned to 154 epiclusters. Viral genomic data was available for 128 cases. Public health authorities identified two large epiclusters (280 and 101 cases each) with a central role on the spread of the disease. Still, the genomic data revealed that each epicluster included two distinct SARS-CoV-2 genetic profiles and thus more than one transmission network. We were able to show that the initial transmission dynamics reconstruction was most likely accurate, but the increasing dimension of the epiclusters and its extension to densely populated settings (healthcare and nursing home settings) triggered the misidentification of links. Genomics was also key to resolve some sporadic cases and misidentified direction of transmission. The epidemiological investigation showed a sensitivity of 70%-86% to detect transmission chains. This study contributes to the understanding of the hurdles and caveats associated with the epidemiological investigation of hundreds of community cases in the context of a massive outbreak caused by a highly transmissible and new respiratory virus.

Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution.

Background and objectives: To understand how organisms evolve, it is fundamental to study how mutations emerge and establish. Here, we estimated the rate of mutation accumulation of SARS-CoV-2 in vitro and investigated the repeatability of its evolution when facing a new cell type but no immune or drug pressures. Methodology: We performed experimental evolution with two strains of SARS-CoV-2, one carrying the originally described spike protein (CoV-2-D) and another carrying the D614G mutation that has spread worldwide (CoV-2-G). After 15 passages in Vero cells and whole genome sequencing, we characterized the spectrum and rate of the emerging mutations and looked for evidences of selection across the genomes of both strains. Results: From the frequencies of the mutations accumulated, and excluding the genes with signals of selection, we estimate a spontaneous mutation rate of 1.3 × 10 -6 ± 0.2 × 10-6 per-base per-infection cycle (mean across both lineages of SARS-CoV-2 ± 2SEM). We further show that mutation accumulation is larger in the CoV-2-D lineage and heterogeneous along the genome, consistent with the action of positive selection on the spike protein, which accumulated five times more mutations than the corresponding genomic average. We also observe the emergence of mutators in the CoV-2-G background, likely linked to mutations in the RNA-dependent RNA polymerase and/or in the error-correcting exonuclease protein. Conclusions and implications: These results provide valuable information on how spontaneous mutations emerge in SARS-CoV-2 and on how selection can shape its genome toward adaptation to new environments. Lay Summary: Each time a virus replicates inside a cell, errors (mutations) occur. Here, via laboratory propagation in cells originally isolated from the kidney epithelium of African green monkeys, we estimated the rate at which the SARS-CoV-2 virus mutates-an important parameter for understanding how it can evolve within and across humans. We also confirm the potential of its Spike protein to adapt to a new environment and report the emergence of mutators-viral populations where mutations occur at a significantly faster rate.

Cryptic Prophages Contribution for Campylobacter jejuni and Campylobacter coli Introgression.

Campylobacter coli and C. jejuni, the causing agents of campylobacteriosis, are described to be undergoing introgression events, i.e., the transference of genetic material between different species, with some isolates sharing almost a quarter of its genome. The participation of phages in introgression events and consequent impact on host ecology and evolution remain elusive. Three distinct prophages, named C. jejuni integrated elements 1, 2, and 4 (CJIE1, CJIE2, and CJIE4), are described in C. jejuni. Here, we identified two unreported prophages, Campylobacter coli integrated elements 1 and 2 (CCIE1 and CCIE2 prophages), which are C. coli homologues of CJIE1 and CJIE2, respectively. No induction was achieved for both prophages. Conversely, induction assays on CJIE1 and CJIE2 point towards the inducibility of these prophages. CCIE2-, CJIE1-, and CJIE4-like prophages were identified in a Campylobacter spp. population of 840 genomes, and phylogenetic analysis revealed clustering in three major groups: CJIE1-CCIE1, CJIE2-CCIE2, and CJIE4, clearly segregating prophages from C. jejuni and C. coli, but not from human- and nonhuman-derived isolates, corroborating the flowing between animals and humans in the agricultural context. Punctual bacteriophage host-jumps were observed in the context of C. jejuni and C. coli, and although random chance cannot be fully discarded, these observations seem to implicate prophages in evolutionary introgression events that are modulating the hybridization of C. jejuni and C. coli species.